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We consider a covariance matrix composed of asymmet- 
ric and free random Levy matrices. We use the results of 
free random variables to derive an algebraic equation for the 
resolvent and solve it to extract the spectral density. For an 
appropriate choice of asymmetry and Levy index {a/2 = 3/4) 
the free eigenvalue spectrum is in remarkable agreement with 
the one obtained from the covariance matrix of the SP500 
financial market. Our results are of interest to a number of 
stochastic systems with power law noise. 

1. The basic concept of independence of commut- 
ing random variables has been generalized by Voiculescu 
to noncommuting ones such as random matrices, using 
the powerful theory of free random variables The 
extension covers formally all stable distributions such as 
the Levy ones. Free random matrices with gaussian fixed 
point, have been applied to many physical problems . 

Recently we have extended the concept of free random 
Levy variables to matrices and suggested that the 
results may be relevant for addressing the issue of noise 
in stochastic systems with power law distributions. The 
latters are encountered in physics (e.g. solar wind data), 
biophysics (e.g. heartbeat data) and finances (e.g. finan- 
cial time series) |^,^. 

A quantitative way of addressing the issue of noise in 
gaussian stochastic systems is through the use of covari- 
ance matrices. In nongaussian systems, this concept is 
usually substituted by the one of covariation |^ or tail 
covariance 0. In this letter, we will show that the con- 
cept of freeness as extended to matrices [||, allows for a 
simple analytical understanding of the eigenvalue struc- 
ture of the covariance matrix, where the underlying noise 
is inherently power law distributed. 

In section 2 we recall some results for free random ma- 
trices, which will be useful in the further part of the pa- 
per, with an emphasis on the Coulomb gas construction. 
In section 3, we define the covariance for random asym- 
metric Levy matrices and derive a closed algebraic equa- 
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tion for the resolvent in the large size limit and for fixed 
asymmetry. We solve it to derive pertinent spectral dis- 
tributions for Levy covariances. In section 4, we show 
that our results may be relevant to the understanding of 
financial spectral correlations, a point of recent interest 
iSl. Our conclusions are in section 5. 



2. The resolvent for an iV x random matrix M is 
generically given by 



G(z) = — 



Tr 



1 



z-M 



(1) 



In the following we consider the unitary ensemble of Her- 
mitean matrices (other ensembles can be studied simi- 
larly) in the large- A^ limit. The averaging is carried with 
the measure g-^Try(M) p^^, standard random ma- 
trices the potential V{M) is usually an even power series 
in M bounded from below. For free random Levy matri- 
ces, the potential is not analytic in M . It is known explic- 
itly in two cases g: V{M) = hi{M'^ + 6^) (Cauchy) and 
V{M) = (1/M-f InM^) (Levy-Smirnov), and implicitly 
as we now discuss. 

The measure given above depends only on eigenvalues 
of the M matrix. The rotational degrees of freedom can 
be integrated out yieldieng the Coulomb-gas representa- 
tion iiofi 



i<j 



where 



V'{X) = 2Re G{\ + iO) 



A,)^ (2) 



(3) 



For free Levy matrices G{z) satisfies an algebraic equa- 
tion in the N ^ oo limit Q {a ^ 1) 



bG"{z) ~ {z ~ a) G{z) + 1 ^ , 



(4) 



in the upper half-plane, and follows by Cauchy reflection 
in the lower half-plane. The parameter b is related to the 
Levy index a, asymmetry and range or scale 7 In 
particular, for symmetric distributions (/? = 0), it reads 



b = —7 e' 



i7r/2 



(5) 
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The marginal case a = 1 is discussed in g. Equation 
can be solved analytically in many cases and nu- 
merically by starting at large z and moving inward, us- 
ing continuity and the asymptotic value of the physical 
branch 



G{z) -> Ijz^hlz 



l+a 



(6) 



The spectral density is p{X) = — ImG(A -I- z0)/7r. For fi- 
nite iV we expect corrections proportional to 1/iVP with 
some positive p, depending on the ensemble. In the fol- 
lowing we shall neglect these corrections. 

To show that this is legitimate and that (j^^ define 
correctly the free random Levy ensemble, we have nu- 
merically calculated the spectral density using a Monte- 
Carlo simulation of the Coulomb gas in the Cauchy and 
Levy-Smirnov case for large, but finite TV. In Fig. Ilia we 
show the results for the Cauchy case and in Fig. ^b we 
show the results for the Levy-Smirnov case. For large 
N , the results asymptote {h/TT)/{\'^ + 6^) (Cauchy) and 
V4A — l/(27rA^) (Levy-Smirnov) shown in sohd fines. 



Other cases can be obtained by solving (|^^ to generate 
the potential, and then Monte-Carlo sampling to gener- 
ate any n-point density. Note that the Cauchy and Levy- 
Smirnov cases are directly accessible by random matrix 
sampling since the measure is explicitly known in these 
two cases (and only these two cases). 
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FIG. 1. Right: Metropolis Monte Carlo simulation of the 

Coulomb gas with a — 1/2 and P — 1 (crosses) compared to 

the FRV result (line). Left: the same for a — 1 and /3 — 



diffusive for < a < 1 and sub-diffusive for 1 < a < 2 
with a = 1 the critical divide. 

The spectrum of C contains important information 
about the character of the correlations. It follows from 
the resolvent W{z) of C. Since (|^) involves a product 
MM'^, the resolvent of C can be related to the resolvent 
of M using Voiculescu's powerful machinery of R and S 
transforms for free random variables [0. Following the 
procedure described in in our case we find that Vl^(z) 
satisfies the transcendental equation 



where w{z) satisfies the multi- valued equation 



z = (1 + w){w + m) , 



(8) 



(9) 



for fixed asymmetry m — T/N, and Levy asymmetry 
parameter /3 = 0. The range of the Levy distribution 
enters the formula only through the prefactor of W{z), 
and the normalization of z. Thus, using the scaling on 
the left-hand side of (^) , one can easily rescale the j — 1 
result to any other value of 7. The phase factor on the 
left hand side of can be deduced from (||). 

The distribution of eigenvalues of the covariance ma- 
trix follows from the discontinuity of w{z), 



p(X) = ^Imw(A + iO). 
ttA 



(10) 



The density of eigenvalues of index a and range 7 satisfies 
the self-affine property ^^(A) = p(A/7^/")/7^/". Its tail- 
behavior is given by 



1 TO 



a/2 , 



(11) 



for A ^ Aniin > and symmetric {(3 — 0) entries. 
The distribution is unimodal and unbounded from above. 
Since the covariance is a square, the tail distribution is 
characterized by an index a/2 since the entries Mu, are 
Levy distributed with index a. 



3. One way of quantifying the fluctuations in a 
stochastic system is through the use of a covariance ma- 
trix C, 



1 



(7) 



where Mu is T x N (asymmetric). The normalization 
in ^ follows from the fact that the "variance" for Levy 
matrices grows like T^/" where a = 2 is the expected dif- 
fusive limit (Brownian). Free Levy ensembles are super- 




2 4 6 2 4 6 

X X 
FIG. 2. Left: Spectral density of FRV with a = 1/2 and 
different asymmetry parameters m. Right: spectral density 
for several indices a, at m = 3.22. 
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For a symmetric Cauchy distribution (a = 1), the 
above analysis yields for the resolvent 



1 / (m+1) =F Y^(m-l)2-4mz 



2(1 + z) 

from which the density of eigenvalues follows 



p{\) = - 



1 x/m 



TT A(A+ 1) 



(A — Aiiiin) ^ , 



(12) 



(13) 



with A,„in = (to — l)^/4m. The spectrum starts at the 
minimum eigenvalue Amin and is unbounded from above. 
In general the minimal eigenvalue is determined by 



A. 



(1 + w*)(to + w*) 



(14) 



where is the negative solution of the second order 
equation 

(a-l)wj-(l-|)(m + l)w, -TO = 0. (15) 

Some typical spectra are shown in Fig. We now illus- 
trate the relevance of these distributions for the covari- 
ance spectrum in a financial market, a point of current 
interest 



4. Recently, it was pointed out that financial co- 
variances are permeated by Gaussian noise in the low- 
lying eigenvalue region with consequences on risk assess- 
ment Here we show that the financial covariance 
spectrum is throughout permeated by free Levy noise. 
For that consider the empirical matrix of relative returns 



Cy = ^Y1 MuMtj with 



T 

{x. 



(16) 



it+l 



Xit) /xio - (M) , 



where xa are the raw returns made of the price xa of 
stock i at time t. Several other definitions are also posi- 
ible, see |l^. The returns are normalized to the ini- 
tial price Xio to insure scale invariance. They carry zero 
mean after subtracting the average relative return (M). 
For the raw prices we will use the daily quotations of 
N — 406 stocks from the SP500 market over the period 
ofT+l = 1309 days from 01.01.1991 tiU 06.03.1996 (ig- 
noring dividends). For these data the matrix asymmetry 
is TO = 1308/406 « 3.22. 

In Fig. H we compare the analytical results following 
from the free Levy covariance to the ones following from 
the SP500 market setting the scale of the free covariance 
to match the one of the market. Fig. ^ left shows the 
distribution of eigenvalues for free Levy matrices with 
index a/2 = 3/4 and asymmetry to = 2, versus the raw 
SP500 data. Fig. |3| right shows the same distribution of 



eigenvalues for free Levy matrices versus the reshuffled 
SP500 data. The reshuffled data follow from a random 
reshuffling of the price series for a fixed stock, destroying 
all inter-stock correlations. Our optimal fit preserves the 
index of the cumulative distribution, albeit for a smaller 
asymmetry (to — 2 < 3.22). It is remarkable that our 
free Levy fit to the spectrum of the covariance suggests 
that the time series of returns is power law distributed 
with index a = 3/2 which is consistent with detailed 
studies of the SP500 financial time series 0. 
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FIG. 3. The spectrum of the historically ordered data 
(left) and the reshufHed data averaged over 20 reshufHings 
(right) for SP500 compared to the FRV result with a = 1.5 
and m — 2. 

Finally, we would like to mention that there exists 
an alternative construction for random matrices with 
spectral Levy disorder, as discussed in [|ll|, where the 
matrix elements are populated from independent , one- 
dimensional stable distributions. As a result, these Ran- 
dom Levy Matrix (RLM) ensembles are not rotational 
invariant, contrary to the Free Levy Matrix (FLM) en- 
sembles discussed here. We note that RLM ensembles 
exhibit interesting correlations between the eigenvectors, 
which are absent (or rather trivial) in the case of the ro- 
tation invariant ensembles such as the GUE, GOE, GSpE 
or FLM ones. For detailed studies of covariances based 
on RLM ensembles we refer to [[12i . 



5. We have shown that the resolvent of covariances 
constructed from free random Levy matrices, obeys an 
algebraic equation for all 1 < a < 2. The ensuing dis- 
tribution of eigenvalues was constructed for a range of 
a's. We have argued that this distribution may be rel- 
evant to a large number of problems where power law 
fluctuations are expected, an example being the SP500 
market. Indeed, we have shown that the spectral density 
of the SP500 financial covariance is in overall agreement 
with a free Levy distribution of index a/2 = 3/4. Free 
Levy noise may be dominant in financial covariances, a 
point of relevance to the central issue of assessing risk in 
finances. 
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